Computation of anomalous scaling exponents of turbulence 
from self-similar instanton dynamics 
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We show that multiscaling properties of developed turbulence in shell models, which lead to 
anomalous scaling exponents in the inertial range, are determined exclusively by instanton dynamics. 
Instantons represent correlated extreme events localized in space-time, whose structure is described 
by self-similar statistics with a single universal scaling exponent. We show that anomalous scaling 
exponents appear due to the process of instanton creation. A simplified model of instanton creation 
is suggested, which adequately describes this anomaly. 



One of the major open problems in hydrodynamic tur- 
bulence is the explanation of small-scale intermittency 
characterized by scaling exponents of structure func- 
tions [ij. Deviation of the scaling exponents from Kol- 
mogorov's theory of 1941 based on dimensional consid- 
erations is called the anomaly. The mechanism responsi- 
ble for this anomaly is still not well understood despite 
of strong theoretical effort, e.g., 0. Anomalous scal- 
ing is explained for a class of linear problems including 
the case of Kraichnan passive scalars using the notion of 
zero modes Q. However, there are essential theoretical 
difficulties in application of this method to the nonlinear 
Navier-Stokes turbulence 

Simplified models may help to reveal the main physi- 
cal mechanism leading to intermittency and its relation 
to nonlinear structure of the system. In this respect, shell 
models (Hj-Q are very successful to reproduce statistical 
properties, which agree qualitatively and quantitatively 
with laboratory experiments for homogeneous isotropic 
turbulence. Also, shell models allow for reliable numeri- 
cal simulation at very high Reynolds numbers. 

In this paper, we show that the anomalous scaling ex- 
ponents in shell models of turbulence can be defined in 
terms of velocity maximums, which split into subgroups 
with self-similar statistics. This splitting allows to recog- 
nize instantons as "elementary particles" of the turbulent 
statistics. Instanons are correlated extreme events local- 
ized in space-time and describing a blowup in inviscid 
limit Q. Instanton structure is described by a single 
universal scaling exponent. We establish the exact rela- 
tion between the scaling anomaly and the process of in- 
stanton creation and derive a general formula for scaling 
exponents for a simplified model of instanton creation. 

We focus on the Sabra shell model of turbulence Q 
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Here u n is the complex shell velocity, which can be un- 
derstood as the Fourier component of the velocity field 
at the shell wavenumber k n = fcoA™, /„ is the forcing 
term restricted to the first shells, v is the viscosity, and 
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FIG. 1. (Color online) Typical turbulent dynamics of the shell 
model in the inertial range. Red squares mark maximums of 
speed amplitudes. Instantons are identified by sequences of 
maximums following in increasing time and shell number. In- 
set shows the universal asymptotic form of the blowup propa- 
gating into unperturbed medium (ideal instanton) for inviscid 
system. 

n = 1, . . . , n max with large n ma x- Traditional choice of 
parameters is A = 2, a = 1, b = c = —0.5, in which case 
the "energy" \u n \ 2 and "helicity" £„R"fc„K| 2 are 
conserved in the inviscid system with no forcing. 

The inviscid shell model possesses the self-similar so- 
lution determined up to phase symmetry by H 

u n (t) = -iu*k~ y ° /(««.(*«. - t)ki- yo ), t < U, (2) 

where the parameters > 0, t* and the function f(u) 
are real and yo = 0.281, see Fig. QJinset). Solution fl2J) 
represents the universal asymptotic form of the blowup 
for the inviscid shell model [H . It describes a stable front 
reaching the smallest scale in finite time and leaving be- 
hind the trail u n oc k~ Vo . Note that no singularity ap- 
pears at t = i* for a particular shell number n, since all 
speeds u n (t*) are finite. 

Fig.Q]shows shell speed amplitudes in turbulent regime 
with a high Reynolds number. In numerical simulation, 
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FIG. 2. (Color online) Structure functions determined in 
terms of shell speeds and speed maximums. Shown are the 
plots of log 2 S p (k n ) (red squares) and log 2 S' p (k n ) (black cir- 
cles) for p — 1,...,6. The plots are shifted vertically to 
demonstrate equivalence of scaling laws in the inertial range. 

we used n max = 34 shells with k = 2~ 6 , the forcing 
fx = 2/ 2 = (5 + 5i)lCT 3 , the viscosity v = 5 x lCT 10 and 
the time interval T = 90000. Numerical data displays 
the well-known multiscaling properties characterized by 
the structure functions 

1 r T 

S p (k n ) e Km - \u n \ p dt oc oc A-" c », (3) 

T^oo 1 J 

where the scaling exponents £ p exhibit nonlinear depen- 
dence on p, Fig. [2J Expression @ is valid in the inertial 
range n/ < n < n v , where rif ~ 4 is the largest number 
of the shell affected by forcing (forcing range n < rif) and 
n v ~ 25 is the smallest shell number affected by viscosity 
(viscous range n > n v ). 

Turbulent dynamics represents a series of bursts, which 
differ in intensity and duration leading to the character- 
istic intermittent behavior. It is argued Q that these 
bursts (called instantons) may have self-similar structure 
described statistically by expression (J2J), where the shape 
f(u) and scaling exponent yo are affected by random 
velocity background. For numerical study of instanton 
statistics, we suggest to identify each instanton with a 
sequence of local maximums v n = max|« n (£)| at times 
t n following in increasing order t no < t no+ \ < ■•• < t ni , 
Fig. [T] No maximums of \u n (t)\ and |« n+1 (i)| are al- 
lowed in the interval t n < t < t n+ \. Each instanton is 
created at some shell number no and either reaches the 
viscous range or annihilates at a shell number rij in the 
inertial range. Using this rule, we group all maximums 
of velocity amplitudes into instantons indexed by a. 

Contribution of self-similar solution ^ to the struc- 



ture function integral in Eq. ([3]) is 




u n \Pdt oc (u.fc"*')*- 1 *- 1 cx vl~ l k-\ (4) 



where v n = max|u„(i)| oc u*k~ yo . The scaling law in 
Eq. (j4]) is independent of the function f(u) and exponent 
yo ■ Assuming self-similar structure of instantons, we can 
use Eq. (j4|) to compute structure functions in Eq. ([3|) as 

S p (fc n )cx T lim (5) 

a 

where the sum is taken over all local maximums (all in- 
stantons a). This suggests the alternative definition of 
structure functions in terms of velocity maximums as 

S' p (k n ) = r lim « K Cp oc \- n ^, (6) 

n a 

with the same scaling exponents £ p . Scaling law ([6]) is in 
full agreement with numerical results, Fig. [2] If the sum 
in Eq. ^ is taken over all stable instantons, which prop- 
agate from the creation shell number no all the way to 
the viscous range, the scaling law ((6]) remains valid with 
almost the same exponents £ p (decreased approximately 
by 0.02). Thus, instantons annihilating in the inertial 
range have very small effect on the scaling law in Eq. ((BJ) 
and will not be included in the sum ^ Q from now on. 
These instantons will be discussed later. 

Statistical properties of the instantons can be under- 
stood by considering the structure functions 

Rp, na {k n ) = lim = V vP, n>n , (7) 

»oo 1 *■ — * 

a(n ) 

where the sum is restricted to stable instantons created 
in the shell no. Log-log plots of the functions R Pl n (k n ) 
are presented in Fig. [3^,. After the vertical translation, 
the curves with equal p and different no collapse to the 
same straight line. The slope of this line is proportional 
to p, see Fig. [3Jd, implying the power law dependence 

R P ,n (kn) = c p , no \( n °- n > y , y = 0.225, (8) 

in inertial range. Linear dependence of the scaling expo- 
nent on p means that multiscaling of structure functions 
in Eqs. Q and ((6]) is not a stastistical property for a par- 
ticular instanton. On the contrary, statistics of instan- 
tons created in the same shell is self-similar with a single 
scaling exponent y. Note that y is close but not equal 
to the scaling exponent yo = 0.28f of the instanton in 
Eq. ([2|) propagating in undisturbed velocity background. 

Self-similarity of instantons can be confirmed by con- 
sidering the probability distribution functions (PDFs) 
determining the probability P ng ^ n (v)dv to sample a lo- 
cal maximum v = max |u n (t)| belonging to the instanton 
created in the shell no- Self-similarity for PDFs implies 

Pn ,n(v)=^ n - no)y Pn ,n (\ (n - no)y v), n > n , (9) 
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FIG. 3. Single-scale statistics of instantons created at the 
same shell, (a) Plots of log 2 Rp,n {kn) for structure functions 
corresponding to instantons created at the shell no- Shown 
are the plots for 8 < no < 17, which collapse after vertical 
translation to a straight line for each p = 1, . . . , 5. (b) Plots 
of (1/p) log 2 R P , 12 (k n ) for p = 1, . . . , 7, which collapse after 
vertical translation to a straight line of slope —y. 
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FIG. 4. Self-similarity of PDFs describing local maximums 
for instantons created at shell no- (a) Renormalized PDFs 
\ in °- n)v P n0!n (\ (n °- n)v v) computed for n = 12 and n = 
12, . . . , 20 collapse to a single curve, (b) Similar plots for 
no = 16 and n = 16, . . . , 22 presented in log-scale to show 
exponential form and good match of the tails. 

which is in full agreement with numerical results, Fig. 3) 
By definitions in Eqs. © and (|7J|, 

n 

Sp(k n ) = k n S Rp—l,n Q (kn) 

n °„ =1 (10) 

- jU-l \^ r . \(n -n)(p-l)y 
~ n / j u p— l,no' < ) 

"o = l 

where we also used Eq. (jHJ). The coefficients c Pi „ are 
expressed from Eq. (fit))) as 

c p,n = k no S p+ i(k m ) — A PJ/ fc„o — (fcno-l)- (11) 

Using the scaling laws (((5]) in Eq. (|11[) yields 

c p ,„ cx A no(1 ~^ +l) (12) 



for fixed p, which is confirmed numerically. We see that, 
due to self-similar structure of instantons, multiscaling is 
attributed to the coefficients c Pj „ in Eq. ((5J depending 
on creation shell number no. 

We showed that scaling properties of turbulence in the 
inertial range are determined exclusively by statistics of 
instantons. In this picture, instantons play the role of 
self-similar "elementary particles" of turbulent dynamics, 
and the anomaly must arise in the process of instanton 
creation. According to Eqs. (I?])-©, characteristic ampli- 
tude of the instanton is given by u n ~ v *\(. n o- n )y with 
the universal exponent y = 0.225, where uq is the cre- 
ation shell number and i>* is the initial amplitude. From 
Eq. Q, we find the characteristic lifetime of the instan- 
ton in shell n as d" 1 /;" 1 . Then the total time occupied 
by instantons in the shell number n is found as 

K X K X = TS' (k n ) cx rfc-<° = T, (13) 

a 

where we used Eq. © with Co = 0. Relation (fl~3"|) was 
also confirmed numerically. It shows that instantons are 
dense in the phase space (n, t) and suggests considering 
dynamics in the inertial range as a "gas" of interact- 
ing instantons. Note that instanton amplitude is inverse 
proportional to its lifetime. Therefore, the observed dy- 
namics with long periods of low activity (intcrmittency) 
is due to weak instantons with long lifetimes. 

Understanding of instanton creation process is a com- 
plicated problem. As we already mentioned, unstable 
instantons, which annihilate in the inertial range, have 
very small influence on scaling exponents. However, nu- 
merical simulation shows that the number of maximums 
corresponding to unstable instantons is not negligible: it 
decreases with n from 30 to 5% of the total number. Un- 
stable instantons, which survive few shells, can be viewed 
as random fluctuations. Unimportance of these fluctua- 
tions for scaling laws suggests that statistics is governed 
primarily by instanton interactions rather than by ran- 
dom background as suggested in [§]. 

Below we propose a very simple model of instanton 
creation leading to multiscaling. Let us consider the PDF 
P n (v) of all instantons in shell n given by 

n 

N n P n {v) = M„ ), (14) 

no = l 

where N n is the total number of instantons (maximums) 
in shell n and M no = N no — A„ _i is the number of 
instantons created in shell n . Our model is based on the 
assumption that instantons are created in self-consistent 
way determined by their distribution in the same shell, 

Pn,n( V ) = ^Pn(^v)- (15) 

This assumption has satisfactory agreement with the nu- 
merical results for the scaling exponent x — 0.21. 



4 



Now we can find the structure functions. Using 
Eqs. ([TBI . ((9]) in the right-hand side of Eq. (|T4]) . we write 

n-l 

N„P n (v) = X x M n P n (X x v) + X y J2 Mn.Pn^.^XH) 

"0 = 1 

= X*M n P n {X x v) + A»JV n _iP n _i(A»«), 

(16) 

where we used Eq. (fT4"|) again in the last equality. Struc- 
ture functions ([6]) can be expressed as 

n r 

S' P (k n )=l^Jv^P n (v)dv. (17) 

Multiplying both sides of Eq. (16]) by (Tk n )~ 1 v p ~ 1 dv and 
integrating yields 

S' p (k n ) = X x ^^S' p {k n ) + A J/(1 ~^~ 1 5p(fc„_i). (18) 

Using the exact result £3 = 1 following from the con- 
stant energy flux condition Q , we substitute p = 3 and 
S' 3 (k n ) oc X~ n in Eq. (0 and obtain 

M n /N n = X 2x {l-X- 2y ). (19) 

Substituting Eq. ([19]) into Eq. ([18]), we obtain 

\y(i-p)-i 

W = I-A^-PHI-A-^) 5 ^"- 1 )' (20) 
This relation implies the scaling law in Eq. ([5]) with 

C P = 1 + y(p - 1) + log A (l - A^ 3 -^(l ~ A~ 2 *0). (21) 
One can see that £3 = 1. The condition £ = requires 
x = (log A (l - A*" 1 ) - log A (l - A- 2 ^))/3 = 0.211, (22) 

which is the value we used above. 

Expression (f2Tj) describes anomalous (nonlinear) de- 
pendence of scaling exponents on p in qualitative agree- 
ment with the numerical data, Fig. [5] The theory of 
instanton creation based on relation (fT5"]) is approximate 
and does not take into account, e.g., dependence on the 
PDFs in neighboring shells. Thus, we did not expect very 
good quantitative agreement. Scaling exponents given by 
formulas ([2T]) . (|22"|) depend on a single parameter y. It is 
interesting that the value y = 0.171 yields the exponents 
£ p , which fully agree with their numerical values for all p, 
Fig- EJ Computations show that y = 0.185 provides very 
accurate values for anomalous exponents of the Glcdzcr- 
Ohkitani-Yamada shell model |f|. Finally, y = 1/3 cor- 
responds to the Kolmogorov scaling £ p = p/3. 

In conclusion, we showed that anomalous statistics of 
turbulence in shell models is determined exclusively by 
instanton dynamics. Instantons represent correlated ex- 
treme events reaching the smallest scale in finite time and 



having self-similar structure described by a single univer- 
sal scaling exponent. We showed that the anomaly in 




_1 1 1 1 1 1 .- 

2 4 6 8 10 

P 



FIG. 5. (Color online) Anomalous part ( p — p/3 of scaling 
exponents. Solid lines show the result of the simplified theory 
of instanton creation with different values of parameter y. 
Error bars represent the result of numerical simulation. 

the inertial range is determined by the process of instan- 
ton creation. We proposed a simplified model of instan- 
ton creation and computed the corresponding anomalous 
scaling exponents. 
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